Spin and orbital angular momentum in gauge theories: 
Nucleon spin structure and multipole radiation revisited 
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We address and solve the long-standing gauge-invariance problem of the nucleon spin structure. 
Explicitly gauge-invariant spin and orbital angular momentum operators of quarks and gluons are 
obtained. This was previously thought to be an impossible task, and opens a more promising 
avenue towards the understanding of the nucleon spin. Our research also justifies the traditional 
use of the canonical, gauge-dependent angular momentum operators of photons and electrons in 
the multipole-radiation analysis and labeling of atomic states; and sheds much light on the related 
energy-momentum problem in gauge theories, especially in connection with the nucleon momentum. 
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The dilemma in separating the nucleon spin — As 
a composite particle, the nucleon naturally obtains its 
spin from the spin and orbital motion of its constituents: 
quarks and gluons. From a theoretical point of view, 
the first task in studying the nucleon spin structure is to 
find out the appropriate operators for the spin and or- 
bital angular momentum of the quark and gluon fields. 
Given these operators, one can then study their matrix 
elements in a polarized nucleon state, and investigate 
how these matrix elements can be related to experimen- 
tal measurements. Disturbingly and surprisingly, after 20 
years of extensive discussions of the nucleon spin struc- 
ture P, 0, H, 0, @| , this first task has never been done, and 
has even largely eluded the attention of the community. 

At first thought, it seems an elementary exercise to de- 
rive the quark and gluon angular momentum operators. 
From the Lagrangian Jz? = -\F^F a ^+-tp(i^D fl -m)ip, 
where B M = d M + igA^ and A M = A*T a (with T a the 
generators of the color SU(3) group), one can promptly 
follow Nother's theorem to write down the canonical ex- 
pression of the conserved QCD angular momentum: 

Jqcd = I d 3 x^-f,ip + J d 3 xip^xX ^tVvj 



d A xE a xA a + d 4 xE ai x x VA ai 



Sq + L„ + Sn 



(1) 



and readily identify the four terms here as the quark 
spin (£ = diag. (<?,<?) and £ x £ = ?£), quark orbital 
angular momentum, gluon spin, and gluon orbital an- 
gular momentum, respectively. However, except for the 
quark spin, all of the other three terms are gauge de- 
pendent, thus have obscure physical meanings. In this 
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regard, it should be noted that the total angular mo- 
mentum is nonetheless gauge invariant (as it must be). 
This can be seen from an alternative, explicitly gauge 
invariant expression 0, @] : 



Jqcd = fd 3 x^~T,il> + J d^i/l'f x - 
d 3 xx x (E a x B a ) 



1^ 



= S q + L' q + J' g . (2) 
This is obtained from Eq.([T]) by adding a surface term, 

d 3 xV -[E a (A a x x)}, (3) 



which vanishes after integration. Since all of the terms 
in Eq.(f2]) are separately gauge invariant, it may seem 
appropriate to identify L' q as the quark orbital angular 
momentum, and J' g as the total gluon angular momen- 
tum. However, a further decomposition of J' g into gauge 
invariant gluon spin and orbital parts is lacking. More- 
over, neither L' q nor J' g obeys the fundamental angular 
momentum algebra, J x J — ij (although J' g does when 
the quark field is absent); hence they cannot be the rel- 
evant rotation generators 0. It has long been assumed 
by the community that the reconciliation of gauge invari- 
ance and the angular momentum algebra is not possible, 
and that gauge invariant, local gluon spin and orbital 
angular momentum operators do not exist Q ■ 

The QED problem revisited — Since QED is also a 
gauge theory, the problems above first emerged there. In 
fact, by simply dropping the color indices, Eqs. |T]) and 
([2]) become exactly the expressions for the electron and 
photon angular momenta, which we denote as: 



Jqed — S e 
= S, 



S~y 



(4) 
(5) 



2 



Eq.© is obtained from Eq.© by adding the same sur- 
face term as in ©, but without the color indices. 

Similarly to the situation in QCD, neither Eq.© nor 
Eq.© is fully satisfactory: On the one hand, the canoni- 
cal angular momentum operators in Eq. ([4]) are what peo- 
ple use familiarly in discussing polarized atomic states 
and radiation, but the gauge dependence of these op- 
erators leads to an uneasy concern about many calcu- 
lations. As one example, the labeling of atomic states, 
which uses eigenvalues of the electron orbital angular mo- 
mentum operator L e — J d 3 xtp^xx A Vi/>, seems gauge de- 
pendent! For another example, the multipole-radiation 
analysis, which employs the multipole- field wavefunction 
constructed with photon spin and orbital angular mo- 
mentum operators, seems again gauge dependent! On 
the other hand, the gauge- invariant operators in Eq.©, 
L' e , and J^,, are not appropriate for constructing angular 
momentum eigenstates (because, as we remarked above, 
they are not angular momentum operators at all), and 
do not separate photon spin from photon orbital an- 
gular momentum. It is stated in common textbooks 
that gauge invariance prohibits the separation of pho- 
ton angular momentum into spin and orbital contribu- 
tions [E IE] > yet both photon spin and orbital angular mo- 
mentum have been measured separately by experiments 

EE El m EE EI El, m . 

Despite the gauge dependence of Eq.©, all QED an- 
gular momentum calculations based on it seem to agree 
well with experiments. It is therefore hard to believe 
that all those discussions, including the whole multipole- 
radiation analysis and labeling of atomic, nuclear, and 
hadronic states, are meaningless. Enlightened by earlier 
clarifications [It], EH , we find that there exists indeed a 
satisfactory and decisive answer for the question of spin 
and orbital angular momentum in QED: 



Jqed = 



1 - f 1 -> 

d 3 X1p / d Z X1p^X X -Dpuretp 



d A xE x A 



phys 



S e + L e + 5L + L . 



d A xE l x x VA 



Here, D pure = V 



xcAp Urei Ap Ure -\- A 



phys 



two parts are defined via: 



V • Aphyg — 0, 



V x A 



pure 



o. 



phys 

(6) 

= A and the 

(7) 
(8) 



These are nothing but the transverse and longitudinal 
components of the vector potential A. The subscripts 
used here are intended to make the physical (vs. pure- 
gauge) content clear, and to prepare for the general- 
ization to QCD. With the boundary condition that A, 
Apure, and A p h ys all vanish at spatial infinity, Eqs. ([7]) 
and © prescribe a unique decomposition of A into A pure 



and Aphys-, and dictate their distinct gauge transforma- 
tion properties: 



A — s. A' —A xUA 

-fipure 7 pure — -^pure T v ji, 



.4 



phys 



A' 



phys 



A 



phys i 



(9) 
(10) 



under a gauge transformation A. Eqs. © and © tell us 
that, in QED, A pure is a pure-gauge field in all gauges, 
and that it transforms in the same manner as does the 
full vector field: A -> A' = A + VA. On the other 
hand, the transverse field Aphys is unaffected by gauge 
transformations, and so can be regarded as the "physical" 
part of A. 

Eq.© is obtained from Eq.© by adding another sur- 
face term, 



d 3 xf ■ [E(A 



■pure 



X X 



)]• 



(ii) 



Now we have all of the elements needed to explain how 
Eq.© gives the correct expressions for the spin and or- 
bital angular momenta of electrons and photons, includ- 
ing their densities. First of all, the total Jqed given by 
Eq.© equals that in Eqs.© and ©, for they merely 
differ by surface terms. Second, the gauge transforma- 
tion properties of A pure and A p h ys show that each density 
term in Eq.© is separately gauge invariant (and hence, 
so is the integrated operator). Third, like the canonical 
L e , the gauge invariant L" satisfies the angular momen- 
tum algebra J x J = ij. This is due to the property 
of A pure in Eq.©. Finally, we note that, in Coulomb 
gauge, V • A = 0, so the longitudinal (pure-gauge) field 
Apure vanishes; thus all quantities in Eq.© coincide with 
their canonical counterparts in Eq.©. This observation 
is of vital importance. It reveals that the gauge invariant 
quantities in Eq. © can all be conveniently computed via 
the canonical operators in Coulomb gauge. This is actu- 
ally what people (implicitly) do in studying atomic and 
electromagnetic angular momenta (such as in multipole 
radiation) , including the recent measurements of t he p ho- 
ton orbital angular momentum [ll|, El, EE Ei , 15, IE]. ^ 
is thus natural that these studies always obtain reason- 
able results. 

Hindsight for QED and solution for QCD — After 
confirming that Eq.© is indeed the correct and satis- 
factory answer for angular momenta in QED, we can ob- 
serve something about it in hindsight: The form of Eq. © 
could have been guessed by reasonable physical consid- 
erations: The photon angular momentum should contain 
only the "physical" part of the gauge field, which should, 
nevertheless, not appear in the expression for the elec- 
tron orbital angular momentum L" e . The latter should 
thus only include the non-physical A pure , so as to cancel 
the equally non-physical phase dependence of the elec- 
tron field, and keep the whole L" gauge invariant. From 
this hindsight for QED, it is natural to expect that the 
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correct, gauge invariant expressions of QCD angular mo- 
menta should be: 

JqCD = [ (fxip^Elp + I (fxip^X X -Bpurelp 



d A xE a x At 



phys 



S a + L" + S" + L 1 ' 



d A xE ai x x VA 



where 



"pure 



igA 



pure 



and 



phys 



(12) 



^pure — ^-pure^ ' The 



essential task remaining now is to properly define the 
pure-gauge field A pure and the "physical" field A phys = 
A phys T a so that they have the desired gauge transforma- 
tion properties, and to prove that the sum of the four 
terms in Eq.([12|) equals that in Eqs.JT]) and ©. This, 
however, turns out to be non-trivial. 

The parallel construction of Eqs.® and ([5]) obviously 
does not work in QCD: For one thing, A pure defined by 
= is not a pure-gauge term in QCD; for 



V x A pure 
another, V-i 



^phys 



and VxA s 



are not invariant 



under the SU(3) gauge transformation: 



C/A M C/ t 



.9 



(13) 



To make . 



^pure 



a pure-gauge term in QCD, we require, 



instead of Eq. ©, that 



1) v A 

"pure ^ *^pure 



V x 



^pure 



^pure 



^pure 



0. (14) 



This provides two independent equations for A pure . We 
still need a third equation that plays the same role as 
Eq.([7]) does in QED, so that A p h ys and A pure have the 
required transformation properties: 



^pure 



k' pure = UA pure tf + ~UVU\ (15) 



phys 



phys 



UA phys U^ 



(16) 



To seek this third equation, we proceed inversely by 
applying these transformations to examine the gauge in- 
variance of each operator in Eq. (|12p . The reason why 
this is possible will be clear shortly below. 

The quark orbital angular momentum L" q provides 
no further constraints. Eqs. (fT4"|) and (fT5j> guarantee its 
gauge invariance, as well as the correct angular momen- 
tum algebra L" q x L' q = iL" q . The gluon spin S" provides 
no further constraints either. Eq. (fT6|) tells us that it is 
gauge invariant. However, the situation for the gluon or- 
bital angular momentum L" g is different: Unlike in QED, 

^phys here is gauge covariant instead of invariant, which 
leads to the gauge transformation of L" : 



E al x x VAX, = 2Tr{E 4 x x VA^ S 



} 



-> 2Ti-{UE l U t xxV(UA phys U ji )} 
= 2Tr{E l f x VA phys } 

+ 2Ti{x x U\VU)(A phys • E - E • A phys )}, (17) 



where E = E a T a . Hence, to make L" g invariant under 
arbitrary gauge transformations, we have to set 



[A phys , E] = Ap h y S • E - E • Aphys = 0. 



(18) 



This is the third equation that we seek. The remaining 
task is to cross-check the consistency of whether or not 
Eqs. fH]) and (Tig)) dictate the transformation properties 
in Eqs.([15l) and p^|) . 

Before making this cross-check, we first make another 



vital check, namely, whether the definitions of . 



^pure 



and 



S^phys by Eqs.fTJJ and (TTg]) ensure that the total an- 
gular momentum in Eq. p2[) is to equal that in Eqs.fl]) 
and @. Since no more tricks are available, the answer 
must be positive or our entire approach will founder. 
A slightly lengthy but straightforward calculation shows 
that the answer is indeed positive: With the definitions in 
Eqs.|U]) and (jTHJ) , Eq.([Tj|) can be obtained from Eq.(T]) 
by adding a surface term similar to (fTlj) for QED: 



d 3 xV-[E a (A; ure xx)}. 



As to the cross-check, we note that A' pure 
given by Eqs. (fT5)l and (fTB)) are solutions of 



and 



5' x A' 

pure pure 
k 'phys^' 



= o, 
= o, 



(19) 

phys 

(20) 
(21) 



where W pure = V-i 



' pure and E' = UEW. The remain- 
ing question is whether Eqs. (f2"0")) and (|2"Tj) have any other 
solution than that given by Eqs. (fT5"|) and (|16[) . This is 
equivalent to asking whether Eqs. (fT4]) and (|T8j) uniquely 
determine the decomposition of A into A pure and A p h ys , 
or, essentially, whether the constraint, [A, E] = 0, fixes 
the gauge completely. This is a tricky question, for, un- 
like in QED, many gauges in QCD suffer from topolog- 
ical complexity such as Gribov copies [19j. Fortunately, 
such complexity does not bother us here: If supplemen- 
tary conditions are needed to restrict the solutions of 
Eqs.dlil) and $ZT$ to that given by Eqs.([T5]) and HI]), 
they can simply be added, without affecting the equiv- 
alence of Eq. (TT2")) with Eqs.fl]) and (O, and without af- 
fecting the gauge invariance of the angular momentum 
operators we constructed; because these properties rely 
only on Eqs.GD, dD, and (TTg]). 

Remarks and discussion — 1) We have noted that for 
QED in Coulomb gauge, V • A = 0, Eq.© coincides with 
Eq.®. Similarly, for QCD, in the gauge [A,E] = (to- 
gether with possible supplementary conditions to com- 
pletely fix the gauge), Eq. (TT2"|) coincides with Eq.([T]). 
Namely, in actual calculations, QCD shares the same 
nice feature as in QED that the gauge- invariant, phys- 
ically meaningful angular momenta can be conveniently 
computed via their canonical, gauge-dependent counter- 
parts in a "physical" gauge in which the pure-gauge com- 
ponent vanishes. From the QCD equation of motion, 
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V • E = ig[A,E] + g^T a iPT a , we see that the gauge 
[A,E] = says essentially that the (gauge-dependent) 
color charge carried by gluons vanishes. So [A, E] = 
has the sense of a "generalized" Coulomb gauge, for it 
leads to an equation of motion V • E a = g^T a ^, similar 
to Gauss' law in QED. 

2) Our construction guarantees that all angular mo- 
mentum operators transform properly under spatial 
translation and rotation. To figure out how they trans- 
form under boost, we need to carry out the canoni- 
cal quantization procedure (preferably in the "physical" 
gauge in which the pure-gauge terms vanish), and com- 
pute the commutators of the angular momentum op- 
erators with the interaction-involving boost generators. 
This non-trivial task will be our next work. 

3) In the literature, there have been various discussions 
about decomposing the Yang-Mills field into several com- 
ponents representing different degrees of freedom, based 



21, 2^ 



mainly on group-theoretical considerations [2C 
It would be very interesting to investigate how these de- 
compositions are related to ours, which is dictated by the 
requirement of a physically meaningful angular momen- 
tum expression. 

4) The so-called gluon polarization AG being measured 
at several facilities Q is related to S g in Eq.(fT]) in the 
temporal gauge in the infinite-momentum frame of the 
nucleon [23j. From our discussion, AG is not the gauge 
invariant gluon spin S' g ' that we construct here. 

5) Beth made a direct measurement of the photon spin 
over 70 years ago • Detection and manipulation of the 
photon orbital angular momentum have also been car- 
ried out recently, and have become a hot topic due to 
their potential application in quantum information pro- 
cessing [ll|, [l2|, HI Q EH EH- These measurements 
can be straightforwardly interpreted with the operators 
in Eq. (J6J) , via its equivalence to Eq.((4]) in the Coulomb 
gauge. This should encourage the investigation of the pic- 
ture of the nucleon spin in terms of the gauge-invariant, 
physically meaningful decomposition in Eq. (|12p , which is 
completely analogous to Eq.© for QED. Experimentally, 
the free-beam-based photon measurements can certainly 
not be extended to gluons directly, and appropriate (even 
ingenious!) methods for measuring L" q , S' g ', and L" g will 
have to be invented; but the clear physical meanings and 
explicit gauge invariance of these quantities guarantee at 
least that there can be pertinent theoretical calculations 
of them, especially in lattice QCD. 

6) From the correct, gauge-invariant angular momen- 
tum expression in Eq.©, we can read out the correct 
electromagnetic momentum density to be E l V A phys , in- 
stead of the renowned Poynting vector, Ex B. The latter 
actually includes a spin current, and can be unambigu- 
ously distinguished from the purely mechanical momen- 



tum E l VA l phys by delicate measurement 0. In QCD, 

E a x B a leads to a picture that gluons carry half of the 
nucleon momentum on the light-cone [251 ]. This picture 
may therefore need to be revised. Similarly to the situ- 
ation for the angular momentum, the momentum opera- 
tors we propose transform properly under spatial trans- 
lation and rotation, and next, we will study how they 
transform under boost by computing their commutators 
with the boost generators via canonical quantization. 
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